Analytical solutions to the spinor Bose-Einstein condensates 
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We propose a method to analytically solve the one-dimensional coupled nonlinear Gross-Pitaevskii 
equations which govern the motion of the spinor Bose-Einstein condensates. Classes of exact station- 
ary solutions for the F — 1 and F = 2 condensates are obtained, from which the soliton complexes 
are naturally produced and classified. The spin polarizations exhibit parallel or kinked spin density 
waves. Our method is also applicable to other coupled nonlinear systems in quantum as well as in 
classical physics. 
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Most Bose-Einstein condensates (BECs) of dilute 
atomic gases realized thus far have internal degrees of 
freedom arising from spins. The spinor BECs have been 
realized in the optical traps which confine the atoms re- 
gardless of their spin hyperfine states. The direction of 
the spin can change dynamically due to collisions between 
the atoms T"!3|. The spinor BECs exhibit a rich variety of 
magnetic phenomena. They give rise to phenomena that 
are not present in the single component BECs, including 
magnetic crystallization, spin textures as well as frac- 
tional vortices 4-6]. Moreover, there exists an interplay 
between superfluidity and magnetism due to the spin- 
gauge symmetry. A ferromagnetic BEC spontaneously 
creates a supercurrent as the spin is locally rotatedj?!, Isj. 

A spinor condensate formed by atoms with spin F = 
1 is described by a macroscopic wave function with 
three components 4' = (-0+1 , ^n, V'-i ■ The mean-field 
Hamiltonian is expressed asjol llOl 



H = dr 



m— — 1 



* r ^ 2 
C2 I 
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where the spin-polarization vector F = ip^^mnipn with 
F (i = X, y, z) the rotational matrices. The coupling con- 
stants Co = (50 + 2(72) /3 and 62 = (.92 — .9o)/3, with gp re- 
lating to the s-wave scattering length of the total spin-F 
channel as gp — inh^ap/M. Utot — IV'iP +"00^ + iV'-iP 
and V{r) is the external potential. The full symmetry of 
the Hamiltonian is S0{3) x U{1). In the ground state, 
the symmetry is spontaneously broken in several differ- 
ent ways, leading to a number of possible phases (ol- 12 1. 
The energies are degenerate for an arbitrary state and 
its globally spin-rotational states ^ — U'^, where U is 
the three-dimensional rotational matrix in the spin space 
which is expressed by the Euler angles as U{a, /3, 7) = 



We are concerned with the quasi-one-dimensional (ID) 
F = 1 and F = 2 BEC in a uniform external potential 
{V{r) — 0). The dynamical motions of the F = 1 spinor 
wavefunctions are governed by idtipm = ^H/dip^, which 
are explicitly written as the coupled nonlinear Gross- 
Pitaevskii equations (GPEs), 
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dl + contot + C2{^^F'^)FUi^„, (2) 
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where cq = Co/2a5^ and C2 = C2/2a^ are the reduced 
coupling constants with a± the transverse width of the 
quasi-lD system. There is a lot of theoretical works that 
numerically solve Eqs.([2]) or the corresponding stationary 
equations isl- isjl . Nevertheless, it is of great interest if 
analytical solutions can be acquired since they provide us 
with further insights into the correlations of the conden- 
sates and help to explore the possible topological struc- 
tures in details. There are some works that utilized the 
similar transformation to obtain the soliton solutions for 
systems with the assumption of time or spatial variations 
in the coupling constant jl^ - fl9l | . Various approximations 
are employed to study the solitons such as bright and 
dark solitons in the F — 1 spinor BECs [20- 23]. How- 
ever, the general exact solutions for the spinor BECs, 
especially for the F = 2 system, are still absent in liter- 
ature. Attempts in this direction are usually frustrated 
by the tremendous complexity of the coupled nonlinear 
GPEs. The challenges are two- fold: one is the nonlinear 
density-density interactions which is associated with the 
Co term, while the other is the spin-exchange couplings 
between the hyperfine states which is associated with the 
C2 term. 

In this Letter we propose a general method which si- 
multaneously decouples the nonlinear density-density in- 
teractions and the spin-spin interactions in the GPEs ([2]) . 
We obtain classes of the exact stationary solutions with- 
out any approximations or the assumptions of time or 
spatial variations of the coupling strengths. The solu- 
tions are either in real or in complex forms and are ex- 
pressed by combinations of the Jacobi elliptical functions. 
We found that the spin-polarizations exhibit cither par- 
allel or kinked spin density- wave configurations. The po- 
lar states (|F| = 0) with non-uniform atomic density are 
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also obtained. The stationary soliton or the sohton train 
complexes are naturally produced, which are identified as 
polarized or vector solitons, and polar or scalar solitons, 
respectively. Furthermore, we obtain new and more gen- 
eral stationary solutions by rotating the states in the spin 
space. All solutions have been numerically verified to be 



stationary by the dynamical evolution with Eqs.(l2]). 

First we deal with the F — 1 System. For convenience 
we choose in the following discussions h = M = 1 as 
the units, cq and C2, or oq and 02, are treated as free 
parameters. The stationary GPEs ([5]) are written as 



M'/'l = i-^dl + (co + C2)(|?AlP + |1^0P) + (co - C2)|V'-l|']l^l + C2V'od, 



(3) 



/^V--! = \-\dl + (co + C2)(|V'-l|' + I^OP) + (co - C2)|V'l|']l^-l + C2V'oV'l*- 



The periodic boundary conditions, 

V'™(1) = V'm(O), C(l) = C(0). (4) 

is adopted in our deduction. In order to seek the an- 
alytical solutions we decouple spin-spin interactions in 
the Eqs.® by assuming V"*™ = ^V'm; S = ±1. It is 
easy to show that S — +1 corresponds to (partially) spin- 
polarized states (|F| ^ 0) whereas S = —1 corresponds to 
spin-unpolarized or polar states (|F| — 0). At the same 
time, the nonlinear density-density interactions between 
different components are decoupled by using the unique 
properties of the Jacobi elliptical functions. There are 
real and complex forms of solutions which are explicitly 
described as follows. 

Real solutions. Consider real solutions of the sn-cn-sn 
form 
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'01(2^) \ / Asn{kx,m) 
ipo{x) = Dcn{kx,m) 
ip-i{x) I \5Asn{kx,m) 



(5) 



where sn, cn, (or dn below) are the Jacobi elliptical func- 
tions and k = AiK{m) with m the modulus (0 < to < 1). 
A and D are the real amplitudes. In the context we al- 
ways take the number of periods j = 2 as examples. One 
has 

l^-ip^lV'iP, IVoP-^'-^lV^ip. (6) 

Substituting ^ into the stationary equations ([3]) , we ob- 
tain three decoupled equations 



llm'ip7n = --;^i'm+lm\ll>ra\^1pra, (to = 0, ±1) (7) 



FIG. 1: (color) (a) The density distribution of the stationary 
solutions ((5]) for the F = 1 condensates with 5 = -1-1. m = 0.5 
and ^ = 600. (c) The spin-rotational solution of ((5|. (b, d) are 
respectively the corresponding soliton solutions with m — >■ 1 
and ^l = 8000. 



ing constants 7™ are defined as 



Mil = M - (co + C2 + C25)D'^ 
Mo M - 2(co + C2 + C25)A^ 
7±i = 2co - (cq + C2 + C2(5)^ 

70 = Cq - 2(co + C2 + C25)^. 



(8) 



where the effective chemical potentials jlm and interact- 



The effective chemical potentials jlm and amplitudes 



3 





X X 



FIG. 2: (color) (a) The density distribution of the stationary 
solutions of cn-sn-cn form for the F = 1 condensates with 5 = 
— 1. m = 0.5, /i = 6000. (c) is the spin-rotational solution. 
(b,d) are respectively the corresponding soliton solutions with 
rrn- 1, /i = 80000. 



A, D can be obtained self-consistently as, 
Ao-ifc'(l-2m2) 

^ 71 2(co+C2)+2c2<5 

70 (co+C2)+C25 

Formula © implies that the effective nonlinear interac- 
tions for tp± must be repulsive (7±i > 0) while for tpo 
attractive (70 < 0) which imposes constraints on the 
coupling strengths cq and 02- Here the value of m is 
arbitrary. 

More general solutions are obtained from solution ([5]) 
by making use of the 50(3) rotational symmetry in the 
spin space ■0m — Umnipn- It is clear that ipm are also 
stationary solution. In the following we fix the three 
Euler angles a = Ttt/IB, (3 = tt/4,j = tt/II as examples. 

Figure 1 display the density distribution of the sn-cn- 
sn solution with 6 = -1-1, which is a (partially) spin- 
polarized stationary states. The values of parameters ao 
and a2 are indicated in the figures. Fig. 1(a) illustrates 
the solution ([5]) and Fig. 1(c) is the result after rotation 
of the solution (O in the spin space, where m = 0.5 and 
/i — 600. We note that ■0-1-1 and ■0_i are staggered af- 
ter the spin-rotation. The soliton solutions are naturally 
produced by taking the limit m — > 1. Fig. 1(b) exhibits 
a typical dark-bright-dark soliton (/i — 8000). The spin- 
rotating solution in Fig. 1(d) reveals more complicated 



FIG. 3: (color online) The phase distribution for the F — 1 
complex solutions, (a) n = 1 with m — 0.29. (b) n = 2 with 
m = 0.92. 



soliton structure. Single soliton solutions can be obtained 
by directly substituting the Jacobi elliptical functions in 
Eq.([5]) with the hyperbolic functions. Other soliton so- 
lutions such as bright-dark-bright complexes can be con- 
structed by the same way. These solitons, being (par- 
tially) spin-polarized (|F| ^ 0), may be properly called 
polarized or vector solitons. 

Figure 2 display the stationary solution of the cn-sn- 
cn form with i5 = — 1, which is a polar state with spin- 
polarization |F| = 0. It should be noted that and 
■0_i remain overlapped after the spin- rotation. In this 
case, the soliton limit exhibits bright-dark-bright com- 
plexes. In comparison to the polarized solitons, these 
spin-unpolarized solitons are called polar or scalar soli- 
tons. 

Complex solutions. We assume ip'^i — ipi (S ^ +1) and 
consider the following cn-sn-cn form of complex solution, 

V'l \ //(x)e'«(-) \ 
■00 = Dsn{kx,m) , (10) 
^-1 J V f{x)e-''^-^ J 



where f{x) = \/ A + Bcn^ {kx,m). A, B and D are real 
constants. One has 

|0_i|^ = |0i^ \^,\^ = D\l + ^)-^\i„\\ (11) 

Substituting pT|) into the coupled GPEs we again 
obtain the decoupled equations ([7]) but with a different 
set of effective chemical potentials and interaction con- 
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FIG. 4: (color) The density distribution of each hyperfine 
state for the = 1 complex solutions. Left column are so- 
lutions from Eg. ([Toll with (a) n = l,m = 0.29, /x = 90, (c) 
n = 2, m = 0.92, fi = 800, and (e) n = 3, m = 0.99, ^J. = 2000, 
respectively. Right column are the corresponding rotating so- 
lutions. 



stants, 

M±i = M - (co + 2c2)^(A + B) 
/io = M-2(co + 2c2)(A + B) 

7±i =2co-(co + 2c2)^ 
70 = Co - 2(co -I- 2c2)^ 

The complex form of solutions to the decoupled Eqs.([7]) 
are self-consistently solved as 

2 , 2 

m k 

A= 3i^(2Ai-(l-2m2)fc2) , (13) 

~ 70 

which leads to 

^ - - 2(c„+C2) (14) 

71 = 270 = 4(co -I- C2). 

One observes that the effective interactions 7^ > for 
this form of complex solution. The phase is 

.W^£^<i«, (15) 

where ai = ±{2fliA'^ - 27iv4^ + P{1 - m'^)AB)i is an 
integral constant. The periodic boundary conditions (|4]) 
require that the amplitudes and phase satisfy 

/(l) = /(0), 0(1) - 0(0) = 2j7r X n, (16) 
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FIG. 5: (color online) The spin-polarization distribution for 
the F = 1 stationary solutions, (a) is the real solution of 
Eg.© with 5 = -1-1. (b)-(d) are the complex solutions of 
Eg. pOl) for n = 1, n = 2, and n = 3, respectively. 



respectively, where n is an integer. The periodic condi- 
tion for the phase is satisfied by adjusting the modulus 
m of the Jacobi elliptical functions. 

Figure 3 display the phase distribution for (a) n = 1 
and (b) n = 2, respectively. They exhibit a four-step 
staircase structure for j = 2. Figure 4 show the density 
distribution for each hyperfine state of the complex sta- 
tionary solution (fTO| and its rotational state for (a, b) 
n = 1, (m = 0.29, /z = 90) (c, d) 71 = 2 (m = 0.92,^ = 
800) and (e, f ) n = 3 (m = 0.99, a* = 2000), respectively. 
In contrast to the real solutions, the complex solutions 
bear mass counter-fiows as well as spin currents, yet the 
density distributions remain stationary. 

Figure 5 show the spin-polarization vector distribution 
F(x) of the above stationary solutions. Fig. 5(a) is for the 
real solution ([S]) which exhibits a parallel polarization 
configuration similar to a spin-density-wave. Fig.5(b-d) 
are for the complex solutions (I10[) with n = 1, n = 2 
and n = 3, respectively. The spin-polarizations exhibit 
kinked structures. 

This method also is applicable to the F = 2 con- 
densates. Figure 6 show the sn-cn-sn-cn-sn form of 
real stationary solutions which correspond to the spin- 
unpolarized states. The details of the results will be 
present elsewhere. 

In summary, we have developed a method to systemat- 
ically solve the one-dimensional coupled nonlinear GPEs 
which govern the motion of the spinor BECs. Several 
classes of analytical stationary solutions for the F = 1 
and F = 2 condensates are obtained. Obviously, we can 
also obtain a lot of other stationary solutions with differ- 
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FIG. 6: (color) The real solutions of the sn-cn-sn-cn-sn form 
for the F = 2 condensates, (a) m = 0.5, /i — 90. (b) m — > 
1, ^ = 6000. (c,d) are the corresponding rotating solutions. 



ent combinations of the Jacobi elliptical functions, which 
will be published in another paper. We emphasize the 
our method is general and exact, without any approxi- 
mations, and can be extended to other nonlinear systems 
such as the coupled nonlinear Klein-Gordon equations, 
the dynamical coupled nonlinear Schrodinger equations, 
the coupled KdV equations, or the circularly symmetric 
coupled sexticoscillators, et al. Works in this direction 
are in progress. 
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